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We renormalize the two-nucleon interaction at leading order (LO) in chiral perturbation theory 
using the scheme proposed by Nogga, Timmermans, and van Kolck — also known as modified Wein- 
berg counting. With this interaction, we calculate the energy per nucleon of symmetric nuclear 
matter in the Brueckner pair approximation and obtain a converged, cutofl-independent result that 
shows saturation, but also substantial underbinding. We find that the renormalized LO interaction is 
characterized by an extraordinarily strong tensor force (from one-pion exchange), which is the major 
cause for the lack of binding. The huge tensor force also leads to the unusually large wound integral 
of 40% in nuclear matter, which implies a very slow convergence of the hole-line or coupled-cluster 
expansion, rendering this interaction impractical for many-body calculations. In view of the unusual 
properties of the renormalized LO interaction and in view of the poor convergence of the nuclear 
many-body problem with this interaction, there is doubt if this interaction and its predictions can 
serve as a reasonable and efficient starting point that is improved by perturbative corrections. 
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I. INTRODUCTION 

The problem of a proper derivation of nuclear forces is as old as nuclear physics itself, namely, almost 80 years. The 
modern view is that, since the nuclear force is a manifestation of strong interactions, any serious derivation has to start 
from quantum chromodynamics (QCD). However, the well-known problem with QCD is that it is non-perturbative in 
the low-energy regime characteristic for nuclear physics. For many years this fact was perceived as the great obstacle 
for a derivation of nuclear forces from QCD — impossible to overcome except by lattice QCD. The effective field theory 
(EFT) concept has shown the way out of this dilemma. One has to realize that the scenario of low-energy QCD 
is characterized by pions and nucleons interacting via a force governed by spontaneously broken approximate chiral 
symmetry. This chiral EFT allows for a systematic low-momentum expansion known as chiral perturbation theory 
(ChPT) [1 . Contributions are analyzed in terms of powers of small momenta over the large scale: [Q/ K^Y , where Q 
is generic for a momentum (nucleon three-momentum or pion four-momentum) or pion mass and sa 1 GeV is the 
chiral symmetry breaking scale. The early applications of ChPT focused on systems like tttt [2] and ttN |3|, where 
the Goldstone-boson character of the pion guarantees that the expansion converges. The past 15 years have also seen 
great progress in applying ChPT to nuclear forces [IH^O]. 

However, there is a difference between the purely-pionic and the one-nucleon sector, on the one hand, and two- 
and multi-nucleon systems, on the other hand. Nuclear physics is characterized by bound states that are nonpertur- 
bative in nature. Weinberg showed 31 that the strong enhancement of the amplitude arises from purely nucleonic 
intermediate states ("infrared enhancement"). He therefore suggested a two-step procedure: In step one, ChPT and 
naive dimensional analysis is used to calculate a "potential" which consists of only irreducible diagrams and, in step 
two, this potential is iterated to all orders by inserting it into a Schrodinger or Lippmann-Schwinger (LS) equation 
to generate the amplitude. 

At leading order (LO), the potential consists of static one-pion exchange (IPE) and two non-derivative contact 
terms. At next-to-leading order (NLO), multi-pion exchange starts which involves divergent loop integrals that need 
to be regularized. An elegant way of doing this is dimensional regularization which (besides the main nonpolynomial 
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result) typically generates polynomial terms with coefficients that are, in part, infinite or scale dependent [S]- One 
reason why so-called contact terms are introduced in the EFT is to absorb all infinities and scale dependencies and 
make sure that the final result is finite and scale independent. This is the renormalization of the perturbatively 
calculated NN amplitude (which, by definition, is the "NN potential"). It is very similar to what is done in the 
ChPT calculations of tttt and TriV scattering, namely, a renormalization order by order, which is the method of choice 
for any EFT. Thus, up to this point, the NN calculation fully meets the standards of an EFT and there are no 
problems. The perturbative NN amplitude can be used to make model independent predictions for peripheral partial 

waves pmng. 

For calculations of the structure of nuclear few and many-body systems, the lower partial waves are the most 
important ones. The fact that in S waves we have large scattering lengths and shallow (quasi) bound states indicates 
that these waves need to be treated nonperturbatively. Following Weinberg's prescription [3] , this is accomplished by 
inserting the potential V into the LS equation: 

T{p\p) = V{p',p) + [ f^V{p',p")^^%—T{p",p), (1) 
J (zTTj-^ p^ — p" +ie 

where Mm denotes the nucleon mass. 

In general, the integral in the LS equation is divergent and needs to be regularized. One way to achieve this is by 
multiplying V with a regulator function, e. g., 

V{p',p) V{p\p) e-(p'/^''" e-^P/'^)'" . (2) 

Typical choices for the cutoff parameter A that appears in the regulator are A « 0.5 GeV ^ 1 GeV [TCI [U]. 

It is pretty obvious that results for the T-matrix may depend sensitively on the regulator and its cutoff parameter. 
This is acceptable if one wishes to build models. For example, the meson models of the past [HI [52] always depended 
sensitively on the choices for the cutoff parameters which were, in fact, welcome fit-parameters for achieving a good 
reproduction of the NN data. However, the EFT approach wishes to be fundamental in nature and not just another 
model. 

In field theories, divergent integrals are not uncommon and methods have been developed for how to deal with 
them. One regulates the integrals and then removes the dependence on the regularization parameters (scales, cutoffs) 
by renormalization. In the end, the theory and its predictions do not depend on cutoffs or renormalization scales. 

So-called renormalizable quantum field theories, like QED, have essentially one set of prescriptions that takes care 
of renormalization through all orders. In contrast, EFTs are renormalized order by order, in which case the number 
of adjustable parameters increases. 

As discussed, the renormalization of perturbative EFT calculations is not a problem. The problem is nonperturbative 
renormalization. This problem typically occurs in nuclear EFT because nuclear physics is characterized by bound 
states which are nonperturbative in nature. 

Weinberg's implicit assumption was that the counterterms introduced to renormalize the perturbatively calculated 
potential, based upon naive dimensional analysis ("Weinberg counting"), are also sufficient to renormalize the non- 
perturbative resummation of the potential in the LS equation. Unfortunately, it has turned out that this assumption 
is not quite correct, as pointed out by Kaplan, Savage, and Wise (KSW) [23], and others. The criticism of the 
Weinberg counting scheme resulted in a flurry of publications on the renormalization of the nonperturbative NN 
problem [24-42, . The literature is too comprehensive to discuss all contributions in detail. Let us just mention some 
of the work that has particular relevance to our present paper. 

If the potential V consists of contact terms only (a.k.a. pion-less theory), then the nonperturbative summation, 
Eq. 0, can be performed analytically, which makes it easier to deal with the renormalization issue. However, 
when pion exchange is included, then Eq. ([I]) can be solved only numerically and the renormalization problem is 
less transparent. Perturbative ladder diagrams of arbitrarily high order, where the rungs of the ladder represent 
a potential made up from irreducible pion exchange, suggest that an inflnite number of counterterms is needed to 
achieve cutoff independence for all the terms of increasing order generated by the iterations. For that reason, Kaplan, 
Savage, and Wise (KSW) [13] proposed to sum the contact interaction to all orders (analytically) and to add pion 
exchange perturbatively up to the given order. Unfortunately, it turned out that the order by order convergence of 
IPE is poor in the ^Si-^Di state [24] . The failure was triggered by the singularity of the IPE tensor force when 
iterated to second order. Therefore, KSW counting is no longer taken into consideration (see, however, Ref. [33] )■ A 
balanced discussion of possible solutions can be found in Ref. [53] . 

Some researchers decided to take a second look at Weinberg's original proposal. A systematic investigation and 
re-analysis of Weinberg counting in leading order has been conducted by Nogga, Timmermans, and van Kolck [31] 
in momentum space, and by Pavon Valderrama and one of the present authors (ERA) at LO and higher orders in 
configuration space |301 1321 133) . A comprehensive discussion of both approaches and their equivalence can be found 
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in Refs. [Mllin]- The LO NN potential consists of IPE plus two nonderivative contact terms that contribute only 
in S waves. By numerical calculations^, Nogga et al. find that the given counterterms renormalize the S waves, i.e., 
the naively expected infinite number of counterterms is not needed. This means that Weinberg power counting does 
actually work in S waves at LO (ignoring the m^r dependence of the contact interaction discussed in Refs. ^51 129| ). 
However, there are problems with a particular class of higher partial waves, namely those in which the tensor force 
from IPE is attractive. The first few cases of this kind of low angular momentum are '^Pa, ^P2, and '^I?2- The leading 
order (nonderivative) counterterms do not contribute in P and higher waves, which is the reason for the problem. But 
the second order contact potential provides counterterms for P waves. Therefore, the promotion of, particularly, the 
^Pq and ^P2 contacts from NLO to LO would fix the problem in P waves. To take care of the problem, a fourth 
order contact needs to be promoted to LO. In this way, one arrives at a scheme of 'modified Weinberg counting' [31] 
for the leading order two-nucleon interaction. 

Once cutoff independence of the on-shell NN T-matrix (and NN phase shifts and observables) has been achieved, 
it is of interest to know if cutoff independent results are also obtained when this interaction is applied in nuclear 
few- and many-body systems. Nogga et al. |31j investigated the lightest such system, namely, the three-nucleon 
bound state, and found cutoff independence of the triton binding energy. It is the purpose of this note to conduct a 
similar investigation in heavier nuclear systems. Since finite nuclei are difficult to calculate, we choose nuclear matter 
(infinitely many nucleons). We will show that the renormalized LO two-nucleon interaction leads to converged results 
for the energy per nucleon in nuclear matter. 

Section II briefly describes and repeats the LO renormalization procedure with modified Weinberg counting intro- 
duced in Ref. |31| . In Sec. Ill, we present the novel point of this paper, namely, the calculation of the energy per 
nucleon in symmetric nuclear matter as a function of density and, in Sec. IV, we compare our results with the work 
by other authors. Conclusions are drawn in Sec. V. 



II. RENORMALIZING THE iViV POTENTIAL IN LEADING ORDER 



In naive dimensional analysis ("Weinberg counting"), the order by order expansion of the chiral NN potential is 
given as: 

Vto = +V^rl"^ (3) 
= V^o + + V^l^ + V^l^ (4) 
Vnnlo = V^i^o + Vi^^ +V^l^ (5) 

l^N3LO = + + Vlt^ + V^t^ + V^t^ (6) 

where the superscript denotes the order v of the low-momentum expansion. LO stands for leading order, NLO for 
next-to-leading order, etc.. Contact potentials carry the subscript "ct" and pion-exchange potentials can be identified 
by an obvious subscript. 

The charge-independent one-pion exchange (IPE) potential reads 

4/^ q^ + mi 

where p ' and p designate the final and initial nucleon momenta in the ccntcr-of-mass system (CMS) and q = p ' — p 
is the momentum transfer; (Ti 2 and Ti^2 are the spin and isospin operators of nucleon 1 and 2; 5^, and denote 
axial-vector coupling constant, the pion decay constant, and the pion mass, respectively. We use /^r = 92.4 MeV and 
gA = 1.29 to correct for the Goldberger-Treiman discrepancy. Since higher order corrections contribute only to mass 
and coupling constant renormalizations and since, on shell, there are no relativistic corrections, the on-shell IPE has 
the form Eq. ([T]) in all orders. 

Here, we will specifically calculate LO neutron-proton {np) scattering and take charge-dependence (isospin violation) 
into accoimt. Thus, the IPE potential reads 

V^7\p',p) = -Vx,{m^.) + (-1)^+1 2Fi^(m,±) , (8) 



For the purposes of the present paper, we conduct the discussion in momentum space. We note however, that all information on the 
necessary number of counterterms can be determined a priori and analytically by inspecting the potential in configuration space at 
short distances (30l l32l [33l [36) . 
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FIG. 1: (color online). Phase shifts and mixing parameters of neutron-proton scattering for total angular momentum J < 2 
and Tiab < 350 MeV. The curves display the LO predictions for cutoff parameter A = 0.5 GeV (black dotted), 1 GeV (blue 
dash-dotted), 5 GeV (green dashed), and 10 GeV (red solid). Note that the dashed and the solid curves are, in general, 
indistinguishable on the scale of the figure. The solid dots and open circles represent the results from the Nijmegen multi- 
energy np phase shift analysis [44j and the VPI/GWU single-energy np analysis SM99 respectively. 



where / denotes the isospin of the two-nucleon system and 

We use TO^o = 134.9766 MeV and m^± = 139.5702 MeV gS]. In the LS equation, Eq. 0, we apply 

= ^MpM^ = 938.9182 MeV, (10) 

2 _ M^TubjTub + 2M„) 
^ " (Afp -t- Af„)2 + 2riabAfp ' ^ ' 

where Mp = 938.2720 MeV and Af„ = 939.5653 MeV are the proton and neutron masses [13], respectively, and Tiab 
is the kinetic energy of the incident neutron in the laboratory system ("Lab. Energy"). The relationship between 
and Tiab is based upon relativistic kinematics. 

Besides the IPE potential, Eq. (|8]), the EFT includes contact terms which represent short-range interactions that 
cannot be resolved at the low-energy scale. Furthermore, the contacts are needed for renormalization. Stating the 
contact potentials in partial-wave decomposition, we have one zero-order (z/ = 0) contact in each S wave: 



v}^\'So) = C.5o, (12) 
(35i) = (13) 



Up to this point, we are still applying Weinberg counting. However, as discussed in the Introduction, higher partial 
waves in which the pion's tensor force is attractive need counter terms to achieve cutoff independence — which leads 
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FIG. 1, continued. 



us to modified Weinberg counting. To be specific, two P waves receive counter terms of second order, 

v}^\'Po) = C.p.p'p, (14) 
V}^\'P2) = C.p.p'p, (15) 

and one D wave needs a fourth-order counter term, 

v}t\'D,) = D.o,p'V. (16) 

For the solution of the LS equation, Eq. ([T]), a regulator function is necessary, for which we choose the one given in 
Eq. ([2]) with n = 2. The regulator depends on the cutoff mass A, which we vary over a wide range from 0.5 GeV to 10 
GeV. In S'-waves, we readjust the contact parameter for each choice of A such that the empirical scattering lengths 
{qs = —23.748 fm for ^Sq and at — 5.4170 fm for '^Si) are reproduced. In those P and D waves which carry a contact 
in modified Weinberg counting, the contact parameter is used to fit — for the various choices of A — the empirical phase 
shift at 50 MeV as given in Ref. [H]. In all cases, we then calculate the phase shifts for all energies below 350 MeV. 

The resulting phase shifts and mixing parameters for total angular momentum J < 2 are shown in Fig. 1. The 
curves refer to A = 0.5 GeV (dotted), 1 GeV (dash-dotted), 5 GeV (dashed), and 10 GeV (solid). The curves for 
A = 5 GeV and 10 GeV are, in general, indistinguishable on the scale of the figure, which demonstrates that cutoff 
independence (nonperturbative renormalization) has been achieved. It is known from the work of Nogga et al. |31j in 
momentum space^ that this is possible. Our results represent an independent confirmation. 



^ Again, the analysis is easier in configuration space |30l 1321 133| : see, in particular. Table II and the thorough convergence analysis of 
phase shifts with total angular momentum J < 5 in Ref. | 33| . 
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TABLE I; Partial-wave contact parameters as a function of the cut-off A for the leading-order NN potential in modified 



Weinberg counting. The parameters are defined in Eqs. (12l-(16l 



Partial-wave contact — Cutoff parameter A in units of GeV — 

parameter 0.5 1.0 5.0 10.0 

Ciso (10* GeV-2) -0.109966 -0.087189 -0.06739623 -0.064460345 

C3s^ (10"* GeV"^) -0.076005 1.349900 -0.02692560 0.021786000 

Cap,, (10"* GeV-"*) 0.840321 -0.1722517 0.001856514 0.000384981 

Csp^ {10^ GeV"'') -0.2316105 -0.0700665 -0.00251447 0.001251038 

Dsq^ (10* GeV"*^) -0.3347880 0.3899800 -0.00020581 -0.00001055 

III. NUCLEAR MATTER 

As discussed in the Introduction, once cutoff independence has been achieved for the two-nucleon system, a good 
question to ask is if cutoff-independent predictions are also obtained in the nuclear many-body problem when this 
renormahzed NN potential is applied. Nogga et al. [31] addressed this question for the three-nucleon system where 
they confirmed the cutoff independence of the triton binding energy at LO. We wish to turn to heavier nuclear systems 
and choose nuclear matter as the representative sample. 

By definition, nuclear matter refers to an infinite uniform system of nucleons interacting via the strong force without 
electromagnetic interactions. This hypothetical system is believed to approximate conditions in the interior of heavy 
nuclei. We shall assume equal neutron and proton densities, that is, we will consider symmetric nuclear matter. This 
many-body system is characterized by its energy per nucleon as a function of the particle density. 

We will use the well-established Brueckner-Bethe-Goldstone method (in short: Brueckner theory) [52 HStilS] to 
calculate the nuclear matter energy. In this theory, a central role is played by the Brueckner G-matrix which is a 
solution of the Bethe-Goldstone integral equation 

G{w) = V-V-^G{w), (17) 

where w denotes the starting energy, i7o the unperturbed Hamiltonian, and the Pauli operator Q projects onto 
unoccupied states. In the pair approximation, the energy per nucleon is given by 

— = — (m|t|TO) + ^-^ {mn\G{w)\mn — 7im) (18) 

m<kF rn,n<kF 

where A denotes the number of nucleons, t the kinetic energy operator, and kp the Fermi momentum, which is related 
to the density p of symmetric nuclear matter by 

P=^A- (19) 



The starting energy is chosen on-shell, i.e. 
with single-particle energy 
and single-particle potential 



w = e(m) + e(n) (20) 
e{m) =t{m) + U{m) (21) 



^ 1 E„<fep("^"|G(w)|mn-nm), m < kp ^^2) 
1 0, m > kp 

also known as the "gap" choice for the single-particle potential, since a gap will obviously occur at the Fermi surface. 
The calculations are conducted in partial- wave decomposition and the Brueckner integral-equation is solved by matrix 
inversion, see Ref. f48] for details. 

The Bethe-Goldstone method was originally devised to handle the short distance hard core in nuclear systems. In 
chiral EFT, an expansion both in l//,r and l/M^ is carried out. For short distances r <C l/wjr pion mass effects 
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FIG. 2: (color online). Energy per nucleon, EjA, in symmetric nuclear matter as a function of the Fermi momentum kp 
applying the LO A'^A^ potentials with the various cutoffs used in the phase-shift calculations of Fig. 1. The curve patterns 
represent the same cutoffs as in Fig. 1. 



can be neglected and on purely dimensional grounds an inverse power short distance singularity should be expected: 
to LO V{r) - l/(/^r3), to NLO V{r) - l/(/>^), to NNLO, V(r) - XI^J^Mnt^) etc. It is natural to ask whether 
the G— matrix result converges in the limit A — > oo. In the appendix we show that this is indeed the case for the gap 



choice Eq. ( 22 ) [48] . The proof rests on the finiteness of the off-shell ii'-matrix [49 . This gives some confidence on 
the stability of the numerics for increasing cut-off values. 

In Fig. 2, we display our results for the energy per nucleon in symmetric nuclear matter as a function of density 
(measured by the Fermi momentum kp) applying the LO potentials with the various cutoffs used in the phase- 
shift calculations of Sec. II. The same curve patterns in Figs. 1 and 2 indicate the same cutoffs. The nuclear matter 
curves for A = 5 GeV (green dashed) and A = 10 GeV (red solid) cannot be distinguished on the scale of Fig. 2, 
demonstrating that cutoff independence of the predictions is achieved; or in other words, the red solid curve represents 
the renormalized result which, as expected, is convergent. This curve shows saturation at a Fermi momentum hp k, 1.0 
fm~^ and an energy per nucleon Ej A = —2.6 MeV. 

Based upon various pieces of circumstantial evidence, it is generally believed that the "empirical" saturation prop- 
erties of symmetric nuclear matter are kp — 1.35 ±0.05 fm~^ and EjA = — 16± 1 MeV [52]. Thus, our renormalized 
LO result shows considerable underbinding. In Ref. [31], a triton energy of —3.6 MeV was found for the converged 
LO result which also deviates considerably from the empirical value of —8.5 MeV. 

The chief reason for this lack of attraction is the fact that the tensor force of the renormalized LO interaction is 
unusually strong, as we will explain now. 

A simple indicator for the strength of the tensor force component contained in a given NN potential is the predicted 
i?-state probability of the deuteron, Pp, because the transition from S to D states can only proceed via the tensor 
force. For the LO interaction at A = 5 GeV and 10 GeV, the Pp) comes out to be 7.2% (it's converged). Conventional 
potentials predict Pp) typically lower, namely, between 4 and 6%; for example, the AV18 [HI], CD-Bonn [SU], and 
N3L0 ^ potentials predict 5.76%, 4.85%, and 4.51%, respectively. Historically, the largest Pd ever predicted by a 
"realistic" NN potential was 7.0% by the Hamada- Johnston potential [S2] of 1962. 

In nuclear matter, the so-called wound integral k is known to depend sensitively on the strength of the tensor 
force [32 US] . The wound integral is defined as 

|2 



K = p j \<j,-^A'dT, (23) 

where 4> denotes the uncorrelated two-nucleon wave function and if) the correlated one, which are related by 

G(j3 = Vtp, (24) 

implying 

V' = 0--^G0. (25) 

Hq-W 
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TABLE II: Partial-wave contributions for S and P waves and total contributions to symmetric nuclear matter at a density 
equivalent to a Fermi momentum kp ~ 1.35 fm~^. Unless denoted otherwise; numbers without parentheses or brackets state 
contributions to the potential energy as obtained from using the Brueckner G matrix; numbers in parentheses are corresponding 
results obtained in Born approximation, i.e., for G = V; finally, figures in brackets state contributions to the wound integral. 



Partial wave 


LO" 


N3LO'' 


CD-Bonn 50, 


AV18[55 




-13.42 


-16.74 


-16.76 


-16.07 




(-5.98) 


(-14.73) 


(-12.64) 


(-2.77) 




[0.056] 


[0.008] 


[0.005] 


[0.017] 




-13.54 


-19.42 


-18.96 


-17.10 


'Si 


(+10.65) 


(-12.51) 


(-8.63) 


(+5.99) 




[0.090, 0.112]^^ [0.017, O.Ol?]'^ [0.004, 0.034]^= [0.015, 0.053]'= 




3.24 


3.90 


3.91 


3.88 


'Pi 


(3.27) 


(4.06) 


(4.24) 


(4.23) 




[0.000] 


[0.001] 


[0.002] 


[0.001] 




-1.01 


-3.14 


-3.08 


-3.15 


'Po 


(-6.30) 


(-3.03) 


(-2.41) 


(-2.48) 




[0.109] 


[0.001] 


[0.002] 


[0.002] 




10.17 


9.68 


9.81 


9.74 


'Pi 


(11.09) 


(10.29) 


(11.74) 


(12.08) 




[0.004] 


[0.003] 


[0.006] 


[0.007] 




-5.37 


-7.27 


-7.05 


-6.96 


'P2 


(-1.49) 


(-6.96) 


(-6.25) 


(-6.09) 




[0.015, 0.015]'* [0.001, O.OOl]'' 


[0.003, 0.001]"* [0.002, 0.002]"* 


Total potential 


-22.44 


-37.02 


-36.35 


-33.96 


energy 


(+9.07) 


(-26.72) 


(-17.92) 


(+6.96) 


Kinetic energy 


22.67 


22.67 


22.67 


22.67 


Total 


+0.23 


-14.35 


-13.67 


-11.29 


energy 


(+31.74) 


(-4.05) 


(+4.76) 


(+29.63) 


Total wound 


[0.405] 


[0.050] 


[0.058] 


[0.101] 



"The renormalized LO NN potential of this work with A = 10 GeV. 

''Quantitative N^LO NN potential regularized by a Gaussian with cutoff parameter A = 0.5 GeV |16| . 
'^The ^Si-^Si and the ^Si-^Di contributions to the wound integral are given. 
''The ^P2-^P2 and the ■^P2-^^2 contributions to the wound integral are given. 



The physical significance of the wound integral is that it measures the probability for exiting two nucleons to states 
above the Fermi surface. This probabiUty is large for "hard" and strong-tensor force potentials. According to 
arguments conveyed by Brandow |53j . a n-hole line diagram is proportional to k""^ and, hence, the convergence of 
the hole-line expansion depends on the size of k, with large k, suggesting slow convergence. 

As shown in the bottom row of Table II, the renormalized LO interaction produces a total k of 40.5%, while the 
corresponding numbers are 10.1%, 5.8%, and 5.0% for AV18, CD-Bonn, and N3L0. The Hamada- Johnston potential 
generated a total k of 21.1%. The partial-wave contributions to k listed in Table II (numbers in square brackets) 
confirm that the strong tensor force of the LO interaction is the main reason for the extraordinarily high k. The 
^Si-^Di transition, which depends entirely on the tensor force, contributes 11.2% to the LO k, while it's 5.3%, 3.4%, 
and 1.7% for AV18, CD-Bonn, and N3L0, respectively. In Fig. 3, we show the ^Si-^Di transition potential of the 
LO interaction with A — 10 GeV and of conventional potentials revealing dramatic differences, particularly, for high 
momenta. An unusual difference occurs also in the ^Pq state where LO generates a contribution to k of 10.9%, whereas 
conventional potentials have at most 0.2%. The tensor operator is known to have a large matrix element in the ^Pq 
state. The ^Pq potentials are included in Fig. 3. 

The fact that a strong tensor force (and a large k) leads to less binding energy in nuclear matter and finite nuclei 
can be understood as follows [22]. For the purpose of discussion, let's approximate the Brueckner G by 

G{w) ^Vc-Vt ^^ Vt, (26) 
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FIG. 3: (color online). Half off-shell NN potentials, V/{2Tvf, for partial -wave states as denoted. The potentials displayed are 
Argonne V18 (black dotted), CD-Bonn (blue dash-dotted), N3L0 (green dashed), and the renormalized chiral LO potential of 
Sec. II for A = 10 GeV (red solid). The on-shell point at p = 153 MeV/c (equivalent to Tiab ~ 50 MeV) is marked by a solid 
dot. 



where Vc denotes the central force and Vr = vtSi2 the tensor force component of a given NN potential (with 5*12 
the usual tensor operator). Note now that all quantitative nuclear potentials arc fit to the same NN data and, 
thus, produce essentially the same on-shell T-matrix or, equivalently, the same on-shell if -matrix, which in the above 
approximation is given by 

K{wf) ^Vc- Wt-^Vt , (27) 

t-Wf 

where Wf is the free (purely kinetic) starting energy and V denotes the principal value. 

A potential with a strong Vt (implying a large, attractive second order tensor term) will have a less attractive 
central force Vc to arrive at the same on-shell if-matrix as compared to a potential with a weak tensor force. Now, 



when we enter nuclear matter and calculate the G-matrix Eq. ( 26 1 , the Pauli operator Q [which is absent in the free- 



space Eq. ( 27 1] and a larger energy denominator (due to the single particle potential in the many-body environment) 



reduces the magnitude of the second term of the G-matrix equation. These two medium effects are know as the 



Pauli and dispersion effects. The larger the attractive second order tensor term in Eq. (26), the larger the reduction 
of the attraction through the medium effects. Therefore, potentials which produce large integral terms in the G- 
matrix equation will predict less attraction in the many-body system. When the central force is very strong ("hard" 
potential), the above mechanism applies also to the iterations of the central term. This happens obviously in the ^Sq 
state where no tensor force is involved, but never-the-less a large k occurs for the LO interaction. This is also part 
of the reason why the ^Si-^Si contribution to the wound is large for LO, namely 9.0% (cf. the very hard LO central 
force seen in the ^Si frame of Fig. 3). 



An idea of the size of the integral term in the Brueckner equation, Eq. (17), is also obtained by comparing the Born 
approximation (i.e., G = V) with the full G. We therefore provide in Table II also the Born approximation results 
(numbers in parentheses) for the various partial-wave contributions. 

As explained in length in Ref . [MJ , arguments similar to the above also apply to Faddeev calculations of the three- 
particle energy. Thus, the substantial underbinding of the triton found in Ref. [31] is most likely also related to the 
huge tensor force of the renormalized LO interaction. 

We note that the discussion of the severely reduced attraction in nuclear matter due to "hard" central potentials 
and large tensor forces applies, of course, only to a calculation conducted in the Brueckner pair approximation. The 
huge wound integral suggests that there will be large three, four, and higher hole-line contributions which may provide 
additional binding. However, the evaluation of multi-hole line contributions is extremely involved and cumbersome. 
Similar arguments apply when the coupled-cluster expansion is used to deal with the nuclear many-body problem |55) . 
Thus, an extraordinarily strong tensor force makes it very difficult to obtain converged results in the many-body 
system, which is one reason why a large tensor-force potential is inconvenient, to say the least. 

The best studied phenomenology of nuclear forces is the one-boson-exchange model. This model includes a p meson, 
which produces a tensor force of opposite sign as compared to the pion (cf. Fig. 3.7 of Ref. [H]). Careful studies have 
shown that the reduction of the pion's tensor force at short range by the p-meson is crucial to arrive at a realistic 
strength for the nuclear tensor force [56] . 

In ChPT theory, contributions from heavy mesons, like the p, are too short-ranged to be dissolved, but instead 
contact terms are added to the theory. The set of contacts which appears at NLO (~ Q^) includes a tensor term that 
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may be perceived as simulating p exchange. Therefore, there is a chance that, at NLO or higher order, the problem 
of the extraordinarily large tensor force encountered at LO will be resolved. 

IV. COMPARISON WITH OTHER WORKS AND THE BROADER PERSPECTIVE 

Our results clearly show that the saturation mechanism in nuclear matter is compatible with a non-perturbative 
renormalization of the venerable IPE potential. However, the binding energy turns out to be rather small and the 
approach looks discouraging from a coupled-cluster expansion point of view. Let us therefore analyze our results in 
the light of other chiral approaches to nuclear matter, based in spirit on the EFT concept, where mainly perturbative 
schemes but also low cut-offs have been employed. We will also provide some perspective on future work. 

If the cut-off parameter takes sufficiently small values, perturbation theory becomes applicable, since in such a case 
high momentum components are suppressed. For the smallest cut-off value presented in Fig. 2, A = 0.5 GeV, we see 
a clear attraction which is strongly dependent, in fact linearly, on the Fermi momentum. Actually, chiral symmetry 
based approaches for nuclear matter have also been pursued in Ref. [57] in a purely perturbative scheme. The main 
source of attraction stems from once iterated IFF, is proportional to the density, ~ k^p, and depends linearly also on 
the cut-off parameter, which needs to be fine-tuned to a value A — 0.4 — 0.5 GeV to achieve saturation. The origin 
of the divergence is related to vacuum amplitudes which to n-th order with the IPE potential yield a contribution 
to the T-matrix, Eq. ([l]), scaling by naive power counting as ~ A""^ (pion mass neglected). The present calculation 
contains iterated IPE with additional counterterms to all orders for any value of the cut-off parameter and, as we 
see, it does not exhibit this very strong k^p dependent attraction. Thus, in Fig. 2 only some residual cut-off effect is 
displayed after low energy NN physics has been fixed. Clearly, the separately large perturbative contributions which 
scale with positive powers of A would not converge without inclusion of counterterms to all orders. 

A NLO G-matrix calculation with a finite cut-off of A = 500 MeV and counterterms was also undertaken in Ref. [5S] 
and saturation was found. As noted in [32], there would be a fundamental problem of removing the cut-off at that 
order, since the deuteron becomes unbound due to the strong 1/r^ repulsive interaction in the triplet channels. More 
recently, a perturbative approach has also been proposed where a power-counting scheme is introduced with extremely 
low cut-off values A ^ for which saturation is achieved (SSI [10] . 

Finally, T^ow-k approaches represent a coarse graining of the interaction in the physically accessible NN elastic 
region for CM momenta fc < A « 400 MeV with the result that all high precision potentials fitting data with 

/DOF 1 including the IPE tail collapse into a unique V|ow-k potential. It has been suggested [Blj that low- 
momentum interactions in general and chiral N'^LO interactions (having their own cut-off A 500 — 700 MeV [TB1I17) ). 
in particular, can be treated perturbatively in nuclear matter calculations, as the corresponding Weinberg eigenvalues 
lie inside the unit circle. Given the universality of the Viow-k approach and the fact that all interactions contain IPE 
while successfully describing the phase shifts in vacuum, there seems to be no fingerprint left of chiral dynamics from 
the two-body sector. When three-body chiral forces determined from few-body data are included, realistic saturation 
properties with more controlled uncertainties are obtained [63] . A trading between two and three body forces is 
observed and the role played by 3N forces as the Vjow-k cut-off is increased from A = 400 MeV to A = 560 MeV 
becomes less important although they still produce saturation at realistic densities. 

Clearly, explicit chiral dynamics is enhanced for larger cut-off values but also the theoretical difficulties increase. 
The spirit of the original proposal of Weinberg's was that the potential could be perturbatively defined according to 
a prescribed counting. This is clearly possible at large impact parameters where neither the strength of the nuclear 
force nor our lack of knowledge of its short distance components prevent us from using a fuzzy but sensibly large 
momentum cut-off value. In a sense, the G-matrix approach is close in spirit to this idea, where the strength of 
the nucleon force in the nuclear medium is characterized by an effective interaction. As we have shown, the LO 
IPE interaction yields cut-off independent results, is entirely parameterized by vacuum properties, and has the nice 
feature of saturation. The problem is to define what is meant by NLO, and in nuclear matter a scheme tightly close 
to the standard perturbative EFT idea is expected to face the same problems already found in the vacuum sector and 
described in the Introduction. 

Several possible and perturbatively motivated schemes for the renormalization of the chiral NN interaction sug- 
gested already in Ref. [31 J include considering higher pion exchanges as perturbations around the LO calculation 
using distorted waves. This proposal was analyzed in Refs. _32i i42j, where it was shown to be feasible and cut-off 
independent, but with little gain from a practical viewpoint: many more counter terms were needed and a worse 
description was obtained in the ^Sq channel and the deuteron. Actually, non-perturbative calculations with singular 
potentials behave non-analytically in the coupling constant, i.e. 1//^ with non-integer a j32j . A RG-based program 
for pion- full theories was advanced in Ref. [27], with further detailed results provided in Ref. [28]. This RG analysis 
yields one counter term in each iS-wave of order Q"^, which must be iterated, and one in each ^Pj and ^Dj wave 
of order Q~^^^, which may be iterated. This is very similar to the power counting suggested by Nogga et al. [5T] . 
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Subleading counterterms occur in -'^S'o at Q° and in ^Si at Q^^'^, while the subleading '^Pj and '^Dj terms are of order 
Q"^!^ . Subleading and higher terms ought to be treated perturbatively. Finally each S wave receives an additional 
term at order . The number of terms at in the RG scheme is larger than in Weinberg counting at NNLO; 
however, except for the subleading D-wavc interactions, it is essentially the same as Weinberg counting at N'^LO 

Toy models where a 1/r^ singular potential is perturbed by l/r^ interactions provide useful insight, but the con- 
sequences for more realistic cases have not yet been worked out [3H]. A follow-up KSW scheme has been pursued 
in Ref. [39j . where the short distance singularity is tamed by the introduction of a Pauli-Villars like pion-mass of 
about the /O-meson mass. An important prerequisite for specific calculations is the perturbative renormalizability 
described in Ref. [42]. However, as we discussed in Ref. [36] the huge change needed from the simple LO ^5*0 phase to 
the real observed one does not suggest any sort of small effect in the large momentum region as non-perturbatively 
renormalized calculations suggest. As already mentioned, this partial wave provides an important contribution to the 
energy in nuclear matter. 

Turning to non-perturbative schemes, it was noted in Ref. [49] that off-shell properties of renormalized chiral poten- 
tials do not look much different from more conventional ones as the singularities are effectively removed. Furthermore, 
we know that the inclusion of A degrees of freedom in chiral potentials [lOj lead to a pattern of better convergence, 
where e.g. the NLO-A and NNLO- A deuteron does exist [61], with an acceptable phenomenological success and 
where much larger and naively more natural momentum cut-off values display a better convergence (see also the 
second Ref. [TT] where A ~ IGeV is taken). Actually, this is a case where the discussion on renormalizability be- 
comes pointless as the renormalized and the natural sized cut-offs are not too far apart, since finite cut-off effects 
are less important the more singular the potential |36j . Indeed, at order v in the chiral counting the potential scales 
as l/(A^+^r^+'^) and then the finite cut-off correction, S\{k), to the renormalized phase shifts, Sao{k), behaves as 
Soo{k) — SA^k) = C'(A~^/^^''/^) [35]. Moreover, unlike the A-less renormalized scheme, the deuteron D-state probabil- 
ity becomes Pd 5.8%, a comparable value to conventional and phenomenologically successful potentials. According 
to our discussion, the corresponding wound integral would be small enough as to guarantee a good converging pattern 
of few body correlations in the nuclear many body problem. Thus, a G-matrix approach applied to chiral interactions 
including A-isobar degrees of freedom has some chance of furnishing the multiple theoretical requirements of renor- 
malizability, power counting for the potential (a la Weinberg), and presumably convergence in the nuclear medium. At 
present it is unclear whether such a scheme will be phenomenologically successful as the potentials do not contain the 
important spin-orbit contributions and it remains to be seen if those can be represented by appropriate counterterms. 

V. CONCLUSIONS 

In this work, we have renormalized the two-nucleon interaction at leading order (LO) in chiral perturbation theory 
using the scheme proposed by Nogga, Timmermans, and van Kolck j31j — also known as modified Weinberg counting. 
With this interaction, we have calculated the energy of symmetric nuclear matter in the Brueckner pair approximation. 
We find that the energy per nucleon as a function of nuclear matter density converges to a cutoff-independent (i.e., 
renormalized) result and shows saturation. The predicted value for the energy per nucleon at saturation shows 
considerable undcrbinding, which is in line with the converged LO prediction for the triton binding energy of Ref. [31) . 
We demonstrate that the LO interaction contains an unusually strong tensor force (from pion-exchange) , which is the 
main reason for the lack of binding in few- and many-body systems. In fact, the tensor force is stronger than that of 
any NN potential ever constructed in the 50-year history of realistic nuclear forces. 

The huge tensor force of the renormalized LO interaction leads to the unusually large wound integral of 40% in 
nuclear matter, which implies a very slow convergence of the hole-line expansion and, similarly, the coupled-cluster 
expansion, rendering this interaction impractical for many-body calculations. 

It is well-known from the meson-theory of nuclear forces that the tensor force produced by the pion needs to be 
damped at short range (or high momenta). In conventional models, this is achieved through ttNN form factors and 
contributions from heavy- meson exchange (particularly, p exchange). ChPT, which does not include heavy mesons, 
provides a contact term of tensor structure at NLO. Thus, a more realistic tensor force may be expected at higher 
orders. 

Several possible schemes for the renormalization of the chiral NN interaction that have recently been given serious 
thought are designed to renormalize the LO interaction non-perturbatively (as done in the present work) and to add 
higher order corrections in perturbation theory. However, in view of the unusual properties of the renormalized LO 
interaction and in view of the poor convergence of the nuclear many-body problem with this interaction, there is 
doubt if this interaction and its predictions can serve as a reasonable and efficient starting point that is improved 
by perturbative corrections. To make the interaction more suitable for many-body calculations, one may consider 
to apply renormalization group methods |64i to construct a soft two-nucleon force (2NF) plus a three-nucleon force 
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(3NF) equivalent to the original LO interaction. Then there will be no problem with the application of this "Viow-k" 
2NF in the many-body system. However, it may turn out that the 3NF (needed for proper equivalence with the 
original LO interaction) is so strong that the issue is just shifted towards a problem with the convergence of the 3NF 
contribution. 



Appendix A: Finiteness of G-matrix for singular interactions 



In this appendix, we show the finiteness of the G-matrix for singular interactions such as the IPE potential for the 
gap choice Eq. ( 22 ) [48' , which implies that the potential energy is suppressed as compared to the kinetic energy in 
the large momentum region. The idea underlying the proof is that short distance physics much below the healing 
distance |46) does not depend on the Fermi momentum. It is convenient to start with the extension of the G-matrix 
to any arbitrary energy z. 



and re-write this equation as 



G{z) = V + V- 



Q 



■G{z), 



G(z)-i = - 



Q 



Z-Hn 



(Al) 



(A2) 



The G-matrix used in the main text, Eq. (17) corresponds to taking z ^ w = e(m) + e{n). We now introduce the 
extended ii'-matrix defined as 



K(z) = V + VV^—K(z) , 

z — t 



(A3) 



where V denotes the principal value and t is the kinetic energy operator. The if-matrix used in the main text, Eq. (271, 
corresponds to using z = Wf = t{m) + t{n). Below threshold, the principal value prescription can be removed as the 
pole is never hit. This equation can likewise be written as 



K{z)- 



= V'^ -V- 



1 



t 



(A4) 



In Ref. [49] the finiteness of the ii'-matrix off-shell for short distance singular interactions was established solely on 
the basis of on-shell renormalization conditions. This means that if K is finite on shell then the off-shell extension 
remains finite as well. Note that the renormalization conditions are basically equivalent to a fixed energy, e.g. zero 
energy. Thus, K will remain finite also below threshold. Taking z — w and subtracting K{w)^^ from G(u')^^ we 
have (principal value where ever necessary) 



G{w) = K{w) + K{w) 



Q-l 



1 



1 



Hq w — t W — Hq 



G{w) 



(A5) 



In this equation, the regulator may be effectively removed from the if-matrix at the operator level even if the energy 
w ^ Wf does not correspond to the scattering energy ^49) . If the interaction is attractive the single particle potential 
U{n) < and w < Wf and the (finite cutoff) K(w) involves states below threshold. As we see, the first integral 
only involves states below the Fermi surface and is thus bound and cut-off independent provided fci? ^ A. The only 
remaining piece where the cut-off enters explicitly is in the last too terms involving (w — Hq)~^ — (w ~ t)^^ , which 
for the gap choice, see Eq. (22) [48], vanishes for states above the Fermi surface, making the integral convergent as 
well. Thus, the finiteness of the G-matrix is reduced to the finiteness of the if-matrix ojf-shell. 
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